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(N ■ I. INTRODUCTION 

D ■ Wormholes are topological handles in space-time linking widely separated regions of a single universe, or "bridges" 
joining two different space-times. Interest in these configurations dates back at least as far as 1916 jlj] with punctuated 
revivals of activity following both the classic work of Einstein and Rosen in 1935 || and the later series of works 
initiated by Wheeler in 1955 ||. More recently, a fresh interest in the topic has been rekindled by the works of Morris 
and Thorne ji| and Morris, Thorne, and Yurtsever These authors constructed and investigated a class of objects 
they referred to as "traversable wormholes" . Their work led to a flurry of activity in wormhole physics |J . 

The central feature of wormhole physics is the fact that traversable wormholes are accompanied by unavoidable 
violations of the null energy condition, i.e., the matter threading the wormhole's throat has to be possessed of "exotic" 
properties. The classical matter does satisfy the usual energy conditions, hence wormholes cannot arise as solutions 
of classical relativity and matter. If they exist, they must belong to the realm of semi-classical or perhaps quantum 
gravity. In the absence of the complete theory of quantum gravity, the semi-classical approach begins to play the most 
important role for examining wormholes. However, there are not so much results concerning quantized fields on the 
wormhole background. Recently the self-consistent wormholes in the semi-classical gravity were studied numerically 
in our works || . As well some arguments in favor of the possibility of existence of self-consistent wormhole solutions 
to the semi-classical Einstein equations have been given in works of Khatsymovsky in Ref . || . 

Note that all the mentioned results were obtained within the framework of various approximations, whereas no 
one up to now has succeeded in exact calculations of vacuum expectation values on the wormhole background. The 
reason for this state of affairs consists in considerable mathematical difficulties which one faces with trying to quantize 
a physical field on the wormhole background. To overcome these difficulties, in this work we will consider a simple 
model of the wormhole space-time: the short-throat flat-space wormhole. The model represents two identical copies 
of Minkowski space with excised from each copy spherical regions, and with boundaries of those regions are to be 
identified. The space-time of this model is everywhere flat except a two-dimensional singular spherical surface. Due 
to this fact it turns out to be possible to construct the complete set of wave modes of the massive scalar field and 
calculate the ground state energy. 

The aim of our work is to calculate the ground state energy of the scalar field on the short-throat flat-space 
wormhole background using the zeta function regularization approach | ^0[ , pT| which was developed in Refs. [[L2]-jl4j. 
In framework of this approach, the ground state energy of scalar field <j) is given by 

where 
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is the zeta function of the corresponding Laplace operator. To make the eigenvalues A 2 n ^ discrete we assume the field 
(j> to be put into a large ball with Dirichlet boundary condition. The A? n % are eigenvalues of the three dimensional 
Laplace operator C 

(A =Af n) (3) 

where R is the curvature scalar. 

The expression ([!]) is divergent in the limit s — > which we are interesting in. For renormalization we subtract 
from (Q) the divergent part of it 

E™ = ]im(E(s)-E d ™( S )), (4) 

where 

E div (s) = lim E(s). (5) 

m — >oo 

By virtue of the heat kernel expansion of zeta function is the asymptotic expansion for large mass, the divergent part 
has the following form 



E 



' IV(s) "2(m) (47r)3/ 2 r( S -i) (6) 



x <j B Q m 4 T(s - 2) + B 1/2 m 3 T(s - -) + B 1 m 2 T{s - 1) + B 3/2 mT(s - ~) + B 2 r(s) }> , 

where B Q are the heat kernel coefficients. 

In this case the renormalized ground state energy (^) obeys the normalization condition 

lim E rcn = 0. (7) 

m — >oo 

The organization of the paper is as follows. In Sec. || we describe a space-time of wormhole in the short-throat 



flat-space approximation and analyze the solution of equation of motion for massive scalar field. In Sec. HI we obtain 



close expression for zeta function and ground state energy and calculate corresponding heat kernel coefficients. We 



analyze also the expression for ground state energy for different radius of throat. In Sec IV we discuss our results. 
The Appendixes A and B contain some technical details of calculations. 

We use units h = c = G = 1 (except Sec [TV] ). The signature of the space-time, the sign of the Riemann and Ricci 
tensors, is the same as in the book by Hawking and Ellis fll5||. 



II. A TRAVERSABLE WORMHOLE: THE SHORT-THROAT FLAT-SPACE APPROXIMATION 

In this section we consider a simple model of a traversable wormhole. Assume that the throat of the wormhole 
is very short, and that curvature in the regions outside the mouth of the wormhole is relatively weak. An idealized 
model of such a wormhole can be constructed in the following manner: Consider two copies of Minkowski space, 
M + and with the spherical coordinates (t, r±, 8±, ip±) [Notice: A4+ and M.- have a common time coordinate 

t. One may interpret this fact as the identification t + «-> excise from each copy the spherical region r± < a, 
where a is a radius of sphere; and then identify the boundaries of those regions: (t, a, 8+, <p+) <-» (t, a, 9-, <£-)■ The 
Riemann tensor for this model is identically zero everywhere except at the wormhole mouths where the identification 
procedure takes place. Generically, there will be an infinitesimally thin layer of exotic matter present at the mouth 
of the wormhole. 

Such an idealized geometry can be described by the following metric 

ds 2 =-dt 2 + dp 2 + r 2 (p) (dO 2 + sin 2 6 dip 2 ), (8) 
where p is a proper radial distance, — oo < p < oo, and the shape function r(p) is 

r(p) = \p\+a. (9) 
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It is easily to see that in two regions 7Z+ : p > and 1Z- : p < separately, one can introduce a new radial coordinate 
r± = ±p + a and rewrite the metric (|^) in the usual spherical coordinates: 

ds 2 = -dt 2 + dr 2 ± + r 2 ± {d6 2 + sin 2 6 dip 2 ), 

This form of the metric explicitly indicates that the regions 1Z+ : p > and TZ- : p < are flat. However, note that 
such the change of coordinates r = \p | + a is not global, because it is ill defined at the throat p — 0. Hence, as was 
expected, the space-time is curved at the wormhole throat. To illustrate this we calculate the scalar curvature R(p) 
in the metric (||): 

R{p) = Sa-^ip). (10) 

Let us now consider a scalar field 4> in the space-time with the metric (||). The equation of motion of the scalar 
field is 

(□ - m 2 - £R)(j) = 0, (11) 

where m is a mass of the scalar field, and £ is an arbitrary coupling with the scalar curvature R. In the metric (|^), a 
general solution to the equation (|TT| ) can be found in the following form: 

cf>(t, p,6,L P ) = e- i " t u(p)Yin (0,9), (12) 

where Yi n (6, ip) are spherical functions, I = 0, 1, 2, . . ., n = 0, ±1, ±2, . . . , ±Z, and a function obeys the radial 
equation 

u" + 2^u' + (u 2 - fcll - m 2 -£R^u = 0, (13) 

where a prime denotes the derivative dj dp. In the flat regions 1Z±, where r(p) = ±p + a, r'(p) = ±1, and R(p) = 0, 
Eq.© reads 

.2 ^2 + 



u" + -—u' +[u 2 -m 2 - / ■ .1 u = 0. (14) 
p±a V (P±a) 2 / 

A general solution of this equation can be written as 

uf[\(p±a)} = A±h i l 1) {\(p±a)}+B±h i l 2) [\{p±a)}, (15) 

where 



A = \/ uS 1 — m 2 , \lo\ > to, 

^[^[z] are spherical Hankel functions, and Af , are arbitrary constants. 

The solutions u^[\(p ± a)} have been obtained in the flat regions 7Z± separately. To find a solution in the whole 
space-time we must impose matching conditions for u i ± [A(p ± a)] at the throat p = 0. The first condition demands 
that the solution has to be continuous at p — 0. This gives 

uf[— Xa] = u^[Xa], 



A^h^i-Xa] + BfhPi-Xa] - A+h^[Xa] - B+hf\Xa] = 0. (16) 



To obtain the second condition we integrate Eq.([13|) within the interval (— e, e) and then go to the limit e — > 0. Taking 
into account the following relations 

r(p) = \p\+a, r'(p) = sign/), r"(p)=28(p), 
lim f f(p)S(p)dp = /(0), 
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we find 



du i 



dx 



du. 



x— — Xa 



dx 



x—Xa 



8 ^ +n i 
[Xa] 



(17) 



Substituting Eq. (|T|) into (jl^) gives 



^ (1) '[Aa] + ^U ; (1) [Aa]^ A+ - (h^'lXa] + ^^[Xa^j B+ = 0, (18) 



Aihf[-Xa]+Bihf y [-Xa]- 

where hj^ [±Aa] = (dh^ [x]/dx) x= ±\ a . 

In addition to two matching conditions (|l6|) and ( |l8| ) we must demand regular behavior of the scalar field at the 
infinity. For this aim, we will consider a "box approximation", i.e., we will assume, in an intermediate stage of 
calculations, that the wormhole space-time has an finite radius R, so that \p\ < R, and we will go, in the end, to the 
limit R — + oo. In the framework of the box approximation, we demand that the scalar field becomes to be equal zero 
at the space-time bounds p = ±R. Taking into account Eq. (|l5|) gives 

uf[-X(R + a)] = 0, u+[X(R + a)] = 0, 



or 



A-h ( t 1) [-X(R + a)}+ B~hP[-X(R + a)] = 0, 



(19) 



A+hf ] [X(R + a)] + B+hf ] [X(R + a)} = 0. (20) 

The four conditions (Il6l), (18), ( |l9| ) and (|2(i| ) obtained represent a homogeneous system of linear algebraic equations 
for four coefficients Af^BfT^ As is known, such a system has a nontrivial solution if and only if the matrix of 
coefficients is degenerate. Hence we get 



h[ 1] [-Xa] h[ 2) [-Xa] -h\ x '[Xa] ~h\ z, [Xa 

h^'[-Xa] hf[-Xa] - h f-y[Xa]-^h!p{Xa] -hf[Xa] - ^h™ 

Ad Act 



Ml 



/i ; (1) [-A(i? + a)] h^i-XiR + a)} 












h^lXiR + a)} 



h&lXiR + a)} 



and 



= 0, 



where the functions ^i[A], ^2 [A] arc defined as follows 

iX 



*}[\] = y y/a(a + R) {\ W [X(R + a)]hf ] [Xa] - [X(R + a)]/i, (1) [Xa] } , 



iX a 



Va(a + i?)(^ (1) [A(i? + a )] (^h^[Xa] + hl 2) ' 



[Xa 



h[ 2) [X(R + a)} ( ^^ (1) [Aa] + h^' [Xa 
\ A& 



(21) 



After some algebra one can show that the determinant in the above formula is factorized, and so Eq. ( |2l| ) c & n be 
reduced to the following two relations: 



(22) 
(23) 

(24) 
(25) 



We introduced additional factors in order to simplify formulas that follow. These factors do not change the relations 
(p2|), (p3|). A significance of Eqs. ( p^ ) and ( p3| ) is that they determine a set of possible values of the wave number A, 
i.e., a spectrum for scalar field modes. Resolving Eq. ( p2| ) and Eq. (23) we can obtain two families, respectively: 



\jfc(a,R,Q, pi = 1,2,3,.. 

P2 = 1,2,3,.. 



(26a) 
(26b) 
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III. GROUND STATE ENERGY AND HEAT KERNEL COEFFICIENTS. 



A ground state energy is given by 



E = \T, EE^+V* +m2 ' (27) 

a=l,2 i=0 p=l 

which is, in fact, the zero point ener gy of massive scalar field. This expression is divergent. In the framework of the 
zeta function regularization method |l(]Jll| , the ground state energy is expressed in terms of the zeta function 

E{ S ) = \^Q c (s-)^, (28) 

where 

/ 1 \ oo oo 1/2-s 

& ( s - 5 ) = EE E( 2/ + !) ( A l +- 2 (m) 



a=l,2 ;=o p=l 



is the zeta function associated with Laplace operator £ — A — m 2 — £R. The parameter /i, with dimension of mass, 
has been introduced in order to have the correct dimension for the energy. For simplicity we represent Eq. ( p8] ) in 
slightly different form 



(30) 



where we introduced the function with dimension energy 



which we shall call the zeta function, too 



The solutions A ; ° (a, -B, of Eqs. (p2|), ( p3| ) cannot be found in closed form. For this reason we use the method 
suggested in Ref. |l2[ , which allows us to express the zeta function in terms of the eigenfunctions. The sum over p 
may be converted into a contour integral in a complex A-planc using the principal of argument, namely 



c H) - £ E E<* + 1> J> <* 2 + - 2 )""* £ "-f w. 



(32) 



where the contour 7 runs counterclockwise and must enclose all solutions of Eqs. (p2|), (p3|). Shifting the contour to 
the imaginary axis, we obtain the following formula for the zeta function 

Cfs- 1 -) =-m 2 *™ £ f> + l) f ft^-^Ah^], (33) 
\ Z S n a=i,aj=o Jm 

where the functions (pi]), ( p5[ ) in the imaginary axis A = ik read 

^[ifc] = I v [k(R + a)]K v [ka] - K u [k{R + a)]I v [ka], (34a) 

*?[**] = U- ^[tfc] + ^{J 1/ [fe(J2 + a)]^[te] -K v [k{R+a)]I' v [ka]}, (34b) 

with 

* = l + ±. 

The expression (^3|) may be simplified in the large box limit R 3> a, which we are interesting in. Let us rewrite 
in the following form 
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* m -U m + aylK . M (l- K$S±$ffi ). (35) 



In the large box limit, the second term in brackets obeys the inequality 



K v [k(R + a)]I v [ka] 2mR 

I u [k{R + a)]K v [ka]^ 1 ' 

and brings exponentially small contribution for ground state energy. 
Therefore, in the limit of large box we have 

Vj[ik] « I„[k(a + R))}K u [ka}, (37a) 

« J„[*(a + i?)] { (e- ^) ^[Ao] + ^[*o]} . (37b) 

In this moment we have to make comment on above formulas. Our approach is valid in the case if the functions WJ™ 
in the imaginary axis do not have zeros in the domain of integration in Eq. (|33|). It gives the restriction for £. The 
function ^j[ik] has no zeros in the imaginary axis, but function ^/f[ik] has simple zeros if £ > -j. Indeed, by using 
recurrent formulas for Bessel's function, let us represent the function \&f[ife] in the following form 

Vftik] = I v [k(a + R)} | U ~ I ~ K v [ka] - ^K^kaA . (38) 

Since the Bessel's functions K v are positive, the expression in brackets may change sign and therefore the function 
^f[ife] may have zeros, if 

f-|-i>°- ^ 



The lowest boundary for £ is 1/4 for I = 0. More precisely, in this case we have 

1 1 

2k ^ a (a + R) 

As far as fc > to, the function ^5>q has simple zero at point k — (4£ — l)/o if 



n = - . , ' „/ " R (l ~ e-^+V) - \ - ^) . (40) 



1 ma 

*>4 + — ■ (41) 

For this reason in the paper we will consider ground state energy for £ < 1/4. In opposite case we have to modify our 
approach. 

Taking into account these formulas we may divide the zeta function, as well as ground state energy (|30|), into two 
parts 

<(s-l)=<r(s-l) + C*(s-l), (42) 

where 



Pa./" 



x A (\n{x v K v [vx\) + In \x v (bK v \vx\ + ™K' v [vx]j }) . (44) 

Here, /3r = m(a + R), f3 a = ma, S = £ — | and v = I + |. 

The first part of the zeta function ( |43"| ) depends only on the size of box with throat R' = R + a and the asymptotic 
structure of the space time. It is exactly twice the expression in the flat Minkowski space time without throat |12] 



G 



calculated for a massive scalar field inside a ball of radius B! with the Dirichlet boundary condition. The factor two is 
very easy explained: we consider scalar field living in the double-sided plane. The second part (|44|) does not depend on 
a boundary and it depends only on the radius of throat a and non minimal coupling £. It contains information about 
the space-time under consideration. The same division of zeta function into two parts has been already observed for 
space-time of the thick cosmic string Jl3| and the space-time of a point-like global monopole [M. Because the first 
part of zeta function (fl3| ) has already been analyzed in great details, we proceed now to consideration the second part 
(II)- 

Both expressions (fl3|) and (jij) and ground state energy ( |30| ) are divergent in the limit s — > which we are interesting 
in. According with renormalization procedure, we have to subtract from regularized expression for ground state energy 
( |30| ) all terms which survive in the limit m —> oo. This procedure corresponds to the subtraction five (three without 
boundary) first terms of the DeWitt-Schwinger expansion p2[-p^| . 

Our goal now is to find in closed form the expansion of zeta function ( |44|) at the point (— |) as power series over 
s (for arbitrary mass). For this reason we use the uniform asymptotic series over inverse index for Bessel functions 
of large index and argument given in Ref . jig] . We subtract from and add to the integrand of Eq. (|4|) the uniform 
expansion of it up to terms proportional to v~ 3 . After subtraction we may tend s — ► 0. Second part, which is the 
uniform expansion of integrand, gives us the pole structure of zeta function. Going in this way (see details in Appendix 
|a|) we obtain the following series for zeta function at the point (—5) 



hit 



= (47r) 3/2^ (B _i) {^ r (- - 2 ) + ~ I ) + b ^ s - !) 

+ b a V2 /3 a T(s - 1) + 6°r(.s)j - ^L- {& a ln/3 2 + fi a [/3 Q ]} , 

fifens - 2) + b^ 2 /3 R T(s - |) + b?P R T(s - 1) 



(4^)3/ 2 (a + i?)r(s- |) 

+ b* /2 p R r(s-h + b?r( s ) 



(45) 



(46) 



Here 



(» 



°3/2 



167r 2 (a + i?) 
= 64tt 3 / 2 



{&«ln/% + Q*[/3«]} 



°2 



8jt 

T 
1 



128^_ 6 4^ + ^-§ 
5 105 



> b 5/2 = 



96£ 4 - 72^ 



n 45 , 93 

21f 2 £ + 

s 16 s 640 



(47) 



(48) 



and 



°2 



b% 

315 ' ' 



-47T 3 / 2 , bf 



5/2 



1 



bM-b?/3* R + b; 



R»2 



16 ^ h« 

y 7 ^ 03/2 



.3/2 



60 



3/2 



(49) 



(50) 



Above expressions (|45|) and ( [46[ ) contain all terms which survive in the limit s — > 0. The details of calculation and 
a closed form for £l a [(3 a ] are outlined in the Appendix The function f2 Q [/3 Q ] tends to a constant for f3 a — ► and 
^ Q [/3 Q ] = -& Q hi/? 2 + V^b% /2 /p a + 0(l//3 2 ) for /9 Q -> 00 (a = a, R). 

Comparing the above expression with that obtained by the Mellin transformation taking the trace of heat kernel 
(in three dimensions), 



2> 
1 



(47r)3/ 2 r( s - i) 

B 3/2 mT(s - -) 



dtt s - 3 / 2 K[t] t ^ Q 



B m 4 T{s - 2) + B 1/2 m 3 T(s - -) + £im 2 r(s - 1) 



(51) 



£ 2 r( s ) 
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we obtain the heat kernel coefficients: 



B 

B\/2 



8tt 



[(a + R) 3 - a 3 ] 



An 3 ^(a + R)' 2 , 



B x = 32tt 



16 , 
a + —TT(a + R), 

o 



B 3/2 = QAtt 3 ' 2 



B 2 = — 



-^3/2 



(52) 



B 



5/2 



8tt 
3 a 

16 7T 3 / 2 

3 a 2 



32 



315 (a + R) ' 



96<f - 72£ 3 + 21£ 2 



45 
16 i 



93 
640 



1 7T 



3/2 



60 (a + i?) 2 



note the difference between Eqs. (|45|), ( |46| ) and Ec 
over inverse mass m — + oo but the formulas (filf), 



Using above scheme we calculated also the coefficient B5/2 which we will need later for analysis. We should like to 

( PH ) . The Eq. (|5l]) is an asymptotic expansion of zeta function 
46] ) are correct for arbitrary mass m and small s — > 0. In fact we 
extracted an asymptotic (for m — > 00) part of zeta function in the form ( |5l| ) and saved a finite part of it. In the 
limit 771 — > 00 the finite part tends to zero and both formulas are in agreement. This is the reason that the function 
n a [p a ] = -b a ln/3 2 + V^5 /2 /% + 0(1/ fa) for (3 a — 00 (a = a, R). 

As far as the space-time under consideration has singular two-dimensional surface S with codimension one, we 
cannot use standard formulas obtained for smooth background and we have to utilize formulas obtained by Gilkey, 
Kirsten and Vassilevich in Ref. jlTj . The heat kernel coefficients (|52|) coincide exactly with that obtained from general 
formulas in three dimensions given in Ref. Jl7[ . We have to take into account that extrinsic curvature tensor of surface 
S is obtained as covariant derivative of the outward unit normal vector N a : 



K af} = V a N . 

For this reason this vector has coordinates N a = (0, ±1, 0, 0) on the spheres p = ±R, and 

2 

trK 



(53) 



(54) 



R + a 

in both cases. In Appendix [X| we found general formulas for arbitrary heat kernel coefficients and traced out them in 
manifest form up to 63. 

To obtain the ground state energy we have to subtract from our expressions (|30|), (|43|), ( ^i| ) all terms which will 
survive in the limit m — > 00. Then we set s = and radius of box R — > 00. Therefore we arrive to the following 
expression 



1 



327r 2 a 



{b a \n(3 2 a + n a [p a }} 



(55) 



A similar general structure for the ground state energy in massless case has been obtained firstly by Blau, Visser and 
Wipf [0 using dimensional consideration only and it was confirmed by detail calculations in Refs. Jl^,Q. 

Using above-mentioned behavior of the f2 a [/3 a ] , the ground state energy tends to zero for large radius of throat 



E ren 

and it is divergent for small radius of throat: 

E ren 



h a 
°5/2 



327T 3 / 2 ma 2 : 



00, 



b a 

ln(ma), a 



167T 2 a 



(56) 



(57) 



The numerical calculations of ground state energy E ren /m (|55|) as a function of f3 a = ma is depicted in Fig. [j] and 
for £ = i and £ = 0, respectively. The details of numerical calculations is analyzed in Appendix [b|. 
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ma 



FIG. 1. The ground state energy E = E ren /m as a function of ma for fixed mass m and £ = |. The energy has minimum 
at point ma « 0.16 with depth J5 min /m pa —0.0025 



E 



ma 



FIG. 2. The ground state energy E = E ren /m as a function of ma for fixed mass m and £ = 0. There is no minimum energy; 
it is always negative. 



IV. DISCUSSION 



We have calculated the ground state energy of the massive scalar field on the short-throat flat-space wormhole 
background (see Eq. ([55])). It can be written down in the formQ 



he 

= -f(/3 a ) 

a 



(58) 



where f3 a — mca/fi, and /(/3 ) is a function of (3 a which has the asymptotic 

/09a) * ^2 ^ Pa, 0a - 0, 



fifia) 



h a 
°5/2 



32^ 3 / 2 /? Q ' 



Pa 



To characterize the behavior of the ground state energy as a function of £ we note that the coefficient b 5 / 2 is positive 
for all values £ and hence the ground state energy tends to —0 while as f3 a — > oo. In the limit (3 a — > 0, the behavior 
of the ground state energy is determined by the sign of &2 (see Eq. (pT\j) and depends on £. For £ < £* « 0.123, the 
62 is negative and the ground state energy tends to minus infinity, otherwise it tends to plus infinity. This difference 
in an asymptotical behavior at f3 a — > results in two qualitatively different pictures describing the behavior of the 
ground state energy. In the first case £ < £*, the ground state energy is monotonically increasing from — oo to and 
has no extremum (see Fig. 2); while in the second case £ > £*, it has a global minimum. For example in Fig. 1 the 



: In this section we use dimensional units G, c, and h. 
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graph of E ren /m versus (3 a is shown for £ = A. It is seen that the ground state energy has the minimum at (3 a ~ 0.16 
with depth E min /m » —0.0025. 

Let us now speculate about the result obtained. Suppose that the quantum massive scalar field plays the role of the 
"exotic" matter maintaining an existence of the short-throat flat-space wormhole in a self-consistent manner. This 
means that the semiclassical Einstein equations have to be satisfied, 

= ^<T^) ran , (59) 
c 4 

where G^ v is the Einstein tensor, and (T^) is the renormalized vacuum expectation values of the stress-energy tensor 
of the scalar field. The total energy in a static space-time is given by 

E = [ ey^^x, (60) 
Jv 

where e — — (T*) ren = — G*c 4 /87rG is energy density, and the integral is calculated over the whole space. In the 
spherically symmetric metric (|8[) we obtain 



Using the relations G\ = 2r" /r + (r' 2 — l)/r 2 and r(p) = \p | + a we can calculate 

B--*£. (62, 

Note that the total energy is negative. 

In the self-consistent case the total energy must coincide with the ground state energy of the scalar field. Equating 



Eqs. (68) and (62) gives 



G a 



a= iJim, (63 ) 



where lp = ^TiG/c 3 is the Planck length. To make further estimations we take into account that in order to be stable 
a quantum system should be in the state with minimum of ground state energy. This requirement can be fulfilled in 
case £ > £*. In particular, for £ = 1/6 the minimum E m i n /mc 2 ~ —0.0025 is achieved at (3 a — raca/fi « 0.16. This 
gives f(P a ) ~ 4 x 10~ 4 , so that 

aw 0.0141 l P , (64) 

and 



m w 11.35 mp, (65) 

where mp = (hc/G) 1 / 2 is the Planck mass. 

Thus, our estimations have revealed that the self-consistent semiclassical wormhole, if exists, should possess the 
throat of sub-Planckian radius, and the quantum scalar field maintaining the wormhole's existence should have super- 
Planckian mass. Of course, it should be noted that our consideration has been restricted by a toy model of the 
short-throat flat-space wormhole, and so one may expect that in more realistic models results would slightly be 
changed. 

Let us emphasize that the result obtained in this work for the wormhole configuration can be generalized. Really, 
the behavior of ground state energy for small ( p7|) and large (]5^) values of the throat's radius a depends only on two 
dimensionless heat kernel coefficients 62 and 65/2, respectively. Instead of the radius a, we could use a typical size of 
system A (throat) and calculate the coefficients 62 and 65/2 on the corresponding background. Now let us consider 
the dimensionless ground state energy E ren /m. Obviously, it will only depend on the dimensionless combination mX, 
and hence the limit of large (small) mass will correspond to limit of the large (small) size of the system. Since for 
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renormalization we have to subtract the fist five terms (up to b 2 ) of expansion for large mass the ground state energy 
in this limit should be proportional to the next non- vanishing term of expansion 



E ren „ 1 1 & 5/2 T( 8 + \) 



2 (4tt) 3 / 2 (mA) 2 r(s - i) 



65/2 (66) 



327r 3 / 2 (mA) 2: 



which coincides with Eq. (|56j). We would like to note that the coefficient b 5 / 2 is non-zero in the limit R — > oo 
for background with singular scalar curvature as it was shown in Ref. fll7|] . For smooth, non-singular geometrical 
characteristics of background it is zero and we have to take into account the next non-vanishing coefficient which is 
63. In this case we have the following expression in the limit mX — > oo: 

The origin of logarithmic term, as well as the behavior for small size of system is following. The structure of poles 
of zeta function does not depend on the parameters of system m and A. The subtraction of the asymptotics for great 
mass brings us the following contribution to the ground state energy 

2 ( JS^H) x { 6o(Am)4r(s - 2) + h ^ {XmfT{s - §> (68) 

+ b 1 (Xm) 2 r( S -l) + b 3/2 (Xm)T( S -^)+b 2 T(s)\ , (69) 

where b a are dimensionless heat kernel coefficients. Taking the limit in this formula we observe that the heat kernel 
coefficients with integer indices will be survived: 

' ' h a (Xm) 4 - hiXm) 2 + b 2 ) hi(Am) 2 . (70) 



327r 2 (Am) \2 
Therefore in the limit A — > one has 

b 2 ln(Am) 

h ~ Tr 2i\ — \ ' V 71 ) 

in agreement with Eq. j57|). 

Therefore the necessary condition that the ground state energy will possess a minimum is following : the coefficients 
b 2 and the next non- vanishing coefficient (b§/ 2 for singular curvature and 63 for non-singular) must be positive. If it 
is so, the discussion above is valid and the self-consistent semi-classical wormhole exist. The radius of throat of stable 
wormhole and the mass of scalar boson in this case depend on the model of wormhole and value of non-conformal 
coupling £. 
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APPENDIX A: 

The uniform asymptotic expansion of the modified Bessel's functions have the form below 

avm - vf 1 . ;.W - ^ £ if. (A!, 

K>A = naz^ y ^i, 1>x] = j*-er> y ^1 

v k=0 y ' fc=0 

where 
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t = 



1 

VT+x 
1 



2 ' 



y/l + x 1 + In 



l + VTTx 2 ' 



= ( J _ 0«fcW + g J (! - 5r)« fc [t]<ft, «o[t] = 1 

« fe +iM = Uk+iM + *(* 2 - 1) 1 5«fcM + *«*[*]} » M*\ = !■ 



(A2) 



Taking into account these formulas in Eq. (H) we obtain power series over s for zeta function. The uniform 
asymptotic expansion (Al) up to z/~™ allows us to take into account terms up to m 3 ~". Because we need for all terms 
which survive in limit m — > oo we use uniform expansion up to n = 3. 

Therefore we have the following expression for the zeta function 



C ( - \ 



2(3 2s COS 7TS 



E 1 

2=0 



02X 1/2- 
2 Pa 



G?X I CE x 



dx 



\n{x v K v [vx]} + \n\x v (sK v [vx\ + ^K' v [vx]^ - ^ {-v) 



-k 



k=-l 



2/? 2s cos tts 



ira 



E< 

i=0 



2\ l/2-» 



§) £ 



fe=-l 



(A3) 



where functions -/V p may be found in closed form for arbitrary index p using simple program in package Mathematica. 
For p > they are polynomial of 3p degree and has the following form 



fe=0 



The first five N p are listed below 



(A4) 



iV = 0, AU = 277, 



N 3 



AS- 1 - 
4 



1 

t H 1 

12 



'-§' 



i 4 - ~t 6 , 



64J3 _ 24(5 2 + 9 5 _ 1^ 
2 64 



t 6 + 



AOS 2 - 25S ■ 



169 
~64~ 



(A5) 



49 329 
Y ~ "64" 



7 179 q 

t 7 H 1 9 . 

576 



We should like to note that the expression ( |A3| ) is identical to original one 
all divergences are contained into the second part. 
Integrating over x with help of integral 



First term is finite in the limit s — > 0; 



dxx I x 2 



2\ 1/2-s 



(1 + X 2 )^/ 2 



r(|- s )r( s + ^) /„ 



2r(? 



p-3+2s 



(A6) 



and taking the limit s — » in the first term we get 



CM -J 



1 



s^O 



16n 2 a 



(47r) 3 / 2 ar(i 



— J2 (-1)%^,]- 



2' k=-l 



(A7) 



where 



+ 2vr)(x) + -JVi 



A/[/3 ] = 32tt^ i> 2 y dx^/a; 2 - ^ ^ ( m K » (vx) + In [^(ra) + °^K' v {vx) 



T N : - 



(A8) 
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a 02^ 1,^ 2(0,1 + s-1) 

A-i = W 2 ns-~ 9 )E r (l + s + 1/2 y 



1=0 



A p = -8^ 2 l3 1 -pf2 

k=0 

2v 



, -2(2k,s + k+- — -). 
T(l + s + 1/2) 2 h 



1=0 



(1 + V 2 /Pl) s \0 a 



(A9) 
(A10) 
(All) 



The first four A p are listed below 



A [AJ =0, 
A 1 [/3 ] = -8tt 



4(5- 



^3 [/?a] 



4tt 3 / 2 

A, 

16?r 



4 /49 329 
35 VT ~ "64" 



2(0,s) + -2(2,s + l) 
6 



1 



:k;(,>- -) 2(i).,. + i)+4^-i)^(2,«+|) + iz(4, a + |) 

169 



(A12) 



64(5 3 - 24(5 2 



19 

64 



£(0,s+l) + - 40(5 - 255 



64 



Z(2,s + 2) 



179 



Z(4 ' S + 3)+ 2520 Z(6 ' S + 4) 



To find the heat kernel coefficients we have to take limit m — » oo in Eq. (A7). The asymptotic expansion of 2(0, q) 
over inverse power of (3 2 was found in Ref . Q : 



z(o, s) = (3 2 a r(s i) + 2^ + S )p- 2l („(~i - 21, -), 

;=o 

where Ch(s, a) is the Hurwitz zeta function 

00 

( H (s,a) = J2( l +«r s > s > !• 



(A13) 



(A14) 



2=0 



Other functions 2(2k, s + k + *tt=) in ( A10) are expressed in terms of Z(0, g) by relation 



Z(2fc,s + fc 



M r(* + ^ + .) p-1 



2 ' ^n!(fc-n)!r(n+2fi + a) v "''" ' 2 
Taking into account above formulas we obtain the following formulas for heat kernel coefficients 

1 



T(s -2 + n 



^2a n ^ p _ 1 (2p - 1), 



p=0 



1 



°n+l/2 



r( s -f + n)^ 



^a„_ p _i(2p), 



(A15) 

(A16) 
(A17) 



where (l,p> 0) 



8?r 



a_ 1 (-l) = — r(a-2), 



a { (-l) = 16tt(-1) 



t C g(-l-2i,§) 
I!(l-2/) 



r( s -i + z), 



^ r(fc + f ) 



(Ais) 
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z-i{p, k) 



; n!(fc-n)! n + ^ 



n=0 



Zl 



k\ r(fc + £ T i + s) p-i 

-r(i + n + H— + s) 



(k-n)lr(n + ^ + s) 



Using these formulas, one obtains the following expressions for heat kernel coefficients 

2 



b% 



b% 



<Stt 
~3 

87T 



K = -— , &i/a = 0, &? = 32^ 



128£ 3 - 64£ 2 



56 68 
!> ~ 105 



^3/2 = 647T 3 / 2 



16 



5/2" 3 



.3/2 



96£ 4 - 72£ 3 + 21£ 2 



45 93 
lfT + 640 



(A19) 



8tt 



4096,. 4096, 4 35584 1088 , 2 848, 144 

£ 5 P H £ 3 £ H 1 

5 ^ 5 ^ 105 ^ 15 ^ 105^ 385 



The coefficient b k / 2 is polynomial of (k — l)-th order over £. The coefficient 62 changes its sign at point w 0.123 
and 65/2 is positive for arbitrary £. 

Our problem now is to take limit s — > in the second part of Eq. (A3). Because the function Z(p,s) with 
p = 2, 4, • • • may be expressed in terms the Z(0, s) only, let us analyze it in details. Let us suppose for a moment 
that f3 a < 1 and represent Z(0, s) as power series over (3 a : 



Z(0, s) = 201 s £ L ^ r (" + «)# n 6r 2n + 2s - 1, - 



(A20) 



The gamma function T(s) has simple poles in points s = 0, — 1, —2, • ■ ■ and the Hurwitz zeta function Ch(s,p) has 
simple pole only in one point s — 1. They have the following expansion near their poles 



r(s - rt) s ^ 



(-1)" 



1 



- + *[n + l] +0(s), Cff(* + l,p)-^o = --*[p] + 0(*) ) 



(A21) 



where ^[x] is the digamma function. 

All divergences of the function Z(0, s) ( |A20 ) in the limit s — > are contained in the first two terms with n = 0, 1. 
The rest of series is finite and we set s = in it. Therefore we obtain the following expression 



Z(0, s) s ^ = 2/3 2s -PlY{s - 1) + -T{s) + -pi [1 - 2 7 - 4 In 2] - - [12&(-1) + In 2] 



24 



12 



(A22) 



2 E 1 f^(2n-l ) ;)+0( S ) 



where Cr( s ) is the Riemann zeta function and 7 is the Euler constant. 

The series in above formula may be anal ytically continued for arbitrary value of [3 a . First of all using series 
representation of Hurwitz zeta function (A14) we represent this series in the following form 



MAO = jr { -^Pl n Cn ( 2n -1 = f> 



-In 1 



n=2 x ' ;=0 

Then, taking into account the integral representation for logarithm 

dt 



(A23) 



In .r 



1 + f 



and the close expression for series below 



00 1 1 
Jo(z 2 ) = X! ^2 +a; 2) = ^2 



vM - — iat I + * 



ix - 2* 



(A24) 
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one has 



J 2 (/3)=/3 ; 



[ dxxl-fy (- 
Jo I V2 



1 



- ixj3 + * - + ix/3 - 2* - 



(A25) 



The function in rhs is analytical in whole complex plane and therefore it gives analytical continuation the series j2(/3 a ) 
for arbitrary value of (3 a . This representation of j2(/3) is preferable for numerical calculations. 
Using the same approach for other Z(p, q) we arrive at the following formulas for ^4fc[/3 ] 



am 



8tt 



7 



1 



} 



Y Pa \ i60 r(s) " ~ PaV{s l) PaT[s 2) f + w - l(/?o) ' 



48 



a 2 = ^0i s | i65 2 /3 a r( s - 1) + A («j - ^ 2 i> + 1) I 

- - y/e {l285 3 - 16<5 2 + 1,5+ ^} r(s) + W3 (/J«), 



+ ^2 (/3a), 



(A26) 



where 



CJ_ 1 (/3 a ) = 87T 

+ /?a 4 
Wl(/3 a ) = -167T 

+ /?a 
^2 (/3a) - -87T 2 



W 3 (/3a) = -y?T 



7 7 1 

-— --C'(-3) + ln2 

160 2 s v ; 240 



2 C'(-l) + lln2 + l 



1 13 2, „ 
-7 h - In 2 

3 ' 36 3 



+ i ln2 +144 



+ 5 



-<(-!)- 3 In 2 



Yi i 




Vl2 7 8 


.lm 2 ) 





1 







2 




16 


K 


Ji,o(/3a) - 2 


Ml 



16 



+ 4 (' 5 -T^)j2(/3a)+g/3S(/3a) 



71 



3360 



, „ 757 71 

In 2 H 1 7 

20160 6720 ' 



+ 5 



1 4 1, n 

—7 h - In 2 

10 ' 5 5 



+ (5 2 [-87 + 8- 161n2] 



(A27) 



+ (5 3 [1281n2 + 64 7 ] 



64<5 3 - 24<5 2 + -5 - — 
2 64 



Jl,l(/3a) + - 



25^+1^ 
64 



+ 35 



49 329 
T ~ "64" 
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j3,l(A») + 42()j4,l(/3o) [> • 



Here we introduced the following notations 

00 



/3 2n Cff I 2n - 3, 



1 



„=3 »(»- jX" -1 ) 



n=l 



(A28) 
(A29) 



The functions j3(/3) and j P , q {(i) with integer p and q are expressed in terms the function jo(/3) only, by relations 



MP) = 


-2/? 4 r 

Jo 


dxx(l — x) 2 


h,m = 






hAP) = 


-2xj (x) 


— x 2 3o( x )\x 


hAP) = 


-3xj {x) 


- ix 2 j' {x) 


kAP) = 


-Axj (x) 


-6x 2 j' {x) 



(A30) 
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For half-integer indexes p and q in Eq. ( A2E ) we use the integral representation for Hurwitz zeta function from Ref. 



and obtain the following formulas 

J h o(0) = -0 

h,M = -0 



Ch(s,v) 



dx 



r(s) Jo 



-(v-l)x xS -l dx 

e x - 1 



(A31) 



o sinh(x) 
dx 



o sinh(a;) 
dx 



sinh(x) 

Substituting the formulas obtained in Eq. ( |A7|) one has 



Ji(20x), 

{iJi{20x) ~20xJ 2 (20x)} , 



20 4 

5Ji(20x) ~ —0xJ 2 (20x) + -0xJ 3 (20x) 



(A32) 



CM* 



(47r) 3 / 2 r(i 



j- (^m 4 r(s - 2) + 6J /2 m 3 r(s - h + b a im 2 T{s - 1) (A33) 
2) I 2 



where 



+ &° /2 mr( s - 1) + ^r(«)} - I JL-n a (/U 



tt Q (/3 a )=^(/3 a ) + ]T (-l) fc ^(/3 a ). 
fc=-i 



(A34) 



In the limit s -> an additional finite contribution appears due to multiplication s In /3 2 and gamma functions in 
brackets (A33). Because the gamma function has simple poles in points 0, — 1, —2, • • •, the heat kernel coefficients 
with integer indexes will give the finite contributions and we arrive at the following expression 



1 



+0 (47r)3/2 a r( S -i) 



+ b a 3/2 a T(s--)+b a 2 T(s) 



- b a o 0*r(s - 2) + b a 1/2 0iT(s --) + bl0iT{s - 1) 



1 



{b a \n0 2 a +n a [0 a }} 



(A35) 



Therefore we divided the zeta function in two parts: the asymptotic singular part of zeta function in standard form 
and finite contribution. 

Using the same approach for C# xt (s — \) one has 



1 



-Ty \b*0%T( S - 2) + b« 2 3 R T(s - |) + fcf^ r(s - 1) 



^0 (47r)3/2( fl + i?)r( s 2 , 

+ &? /2 /^r( s - i) + 6fr( s )} - 167r2( * + fl) {b R \n0% + n*[0 R }} . 



(A36) 



By virtue of the fact that in the limit m — > 00 the above formulas must give us the asymptotic expansion (|5l]) , the 
function fl[0] has the following behavior [0 a — + 00 , (a = a, i?)) 



= -b<* \n0l £ bt l2 0t k n* + § - 2) 

00 , 
= -6 a ln/3 2 +^6« /2 /3^ fe r(--2) 

fc=5 

= -6 a ln/3 2 +v¥^ /2 /3- 1 +0(/3- 2 ). 



(A37) 
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APPENDIX B: 



The m ain problem for numerical calculation the ground state energy is the term given by Eq. (A.8). The 

series in (A8) is low convergent. To calculate this expression let us, first of all, represent it in the following form 



— A 00 

A f = ^ = J2 a »> where 



2=0 



dx 



0/u 

+ 2vr)(x) + -Ni 
v 



v 2 dx 



\nK v {ux) + In 5K v (vx) + —K v (vx) 



(Bl) 
(B2) 



and divide the series in two parts 



N 



A, 



E a »- 

l=N+l 



(B3) 



;=n 



The first sum we calculated numerically. The calculations become lighter due to fact that the Bessel's functions of 
second kind with half-integer indexes are polynomial with simple exponent factor ]l6[ | . In the second sum we use the 
uniform expansion of the integrand over inverse power of index v. Since we have already subtracted first three terms 
Ni, N2 and N3, the uniform expansion of integrand will start from and we obtain the following expression 



l=N+l 



l =N+ l JP/v V V UX P =A 

where N p [t] is the polynomial of 3p degree: 



(B4) 



(B5) 



k=0 



The coefficients a Py k for p = 1, 2, 3 may be singled out from Eq. (A5). 

Then one takes derivatives and integrals in Eq. (B4) and changes the sums over p and I. After this we arrive at 
the following formula 



where 



E °i = J2 n p 

l=N+l p=4 



'.(-lf-pj2a P , k ^ 7V j^h[p,p + 2k,(3,N] : 

01\ -(9-l)/2 



h\p,qJ,N}= J2 v 2 ~ P [ 1 + 

l=N+l 



The function h may be found in closed form for integer p and q. 
The above function h may be estimated by 

h\p,q,(3,N}& (n + 



(B6) 

(B7) 
(B8) 



2-p 



(B9) 



and the series ( |B6[ ) is fast convergent for great N. We use N = 10 and in order to work with precision 1CP 10 it is 
enough to take expansion up to p = 8 (five terms). In fact, this procedure converts the low convergent series to fast 
convergent series over N 2 ~ p . 
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To illustrate above approach we reproduce in Fig. [3] the three steps of calculation of Af. (i) the zero term (thick 
curve), (ii) the contribution of first eleven terms up to I — 10 (middle thickness) and (iii) exact curve (up to p = 8 
in uniform expansion, thin curve). Ten terms I = 1 -j- 10 gives us the correction of 36% for zero term, and the series 
from I = 11 to oo gives us 4% correction. 
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